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1—5 It is long speculated that long duration gamma-ray bursts (GRBs) originate from a 

^ relativistic jet emerging from a collapsing massive star progenitor. 1 4 Although asso- 

^ ciations of core-collapsing supernovae with long GRB afterglows have been identified 

\q in a number of systems, 5 8 including the latest X-ray flash (XRF) 060218/SN 2006aj 

connection detected by Swift, 9 15 direct evidence of a relativistic jet emerging from a 

o : 

^ collapsing star is still lacking. Here we report the detection of a thermal emission com- 

^ ponent (high-T component) accompanying the prompt X-ray emission of XRF 060218, 

Oh; 

q with temperature kT H = 1.2lto 24 keV and effective blackbody radius R H ~ 5 x 10 9 cm. 
This high-T component co-exists with another low-T thermal component as reported 
by Campana et al. 9 for at least 2700 seconds, but evolves independently with respect 

•i-H . 

^ to the low-T component by tracing the lightcurve of the non-thermal component. We 

b ; 

identify this high-T thermal component as the emission of a hot cocoon surrounding 
the relativistic jet, as expected from the theoretical models. 16 18 ' 419 

XRF 060218 was detected with the Swift/Burst Alert Telescope (BAT) on 2006 February 18.149 UT. 
Swift slewed autonomously to the burst and the X-ray telescope (XRT) began to collect data 159 s 
after the burst trigger. We analyze the XRT data in the first orbit (t = 159 — 2740 seconds since 
the burst trigger). We reduce the data with the Swift/XRT tools and process the spectrum with 
the Xspec software package. We first fit the time-integrated spectrum with a combined model of a 
single-temperature (kT) black body model (BB model) normalized to its apparent area at a given 
distance and a single power law (dN(E)/dE oc E~ r ; PL model), incorporating the neutral hydrogen 
absorbtion in both our Galaxy and the GRB host galaxy. The observed effective" black body radius" 
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(R) could be inferred from the normalization term (K) of this model, R B b = K 1 / 2 x -D 10kpc km, where 
-Diokpc is the source distance in units of 10 kpc. The N H of our Galaxy 20 for this burst is ~ 1.1 x 10 21 
cm -2 . We find that the fitting results depend on the initial guess value of kT. If kT is assigned an 
initial guess value less than 1 keV, one can get the best fit results with kT ~ 0.12 keV, the thermal 
component identified by Campana et al.. 9 However, we also find marginally acceptable fitting results 
with a high-T component with kT ~ 1.3 keV, if the initial guess value of kT is assigned greater than 
1 keV. We replace the single power law component with a broken power law model (BKNPL model) 
in our fitting, and find that the high-T component stands out and the two photon indices of the 
broken power law are I\ ~ —1 (fixed) and T 2 ~ 2.08, respectively, with a break energy at E b ~ 0.80 
keV. The photon spectrum with E < E b roughly agrees with a black body spectrum, being consistent 
with the low-T component. This likely indicates that there exists two thermal components in the 
XRT data, with kT ~ 0.1 keV and kT ~ 1.3 keV, respectively. We then use a combined model 
with two thermal and one power law components (BB+BB+PL model) to fit the spectrum. The 
comparison of the fitting results with the BB+PL model and the BB+BKNPL model is shown in 
Table 1. It is found that besides the previous reported thermal component with kT\ = O.llGt ,; ) ) 2 
(low-T component), 9 a new thermal component (high-T component) with kT 2 = 1.21^^1 keV is 
indeed present in the XRT data. Figure 1 shows the accumulated overall X-ray spectrum during the 
first orbit (159-2740 seconds since the burst trigger), with the contributions from the two thermal 
components marked. We derive the effective blackbody radii of the two thermal components as 
R L = 3.2 x 10 12 cm and R H = 5 x 10 9 cm, respectively, suggesting that the two thermal components 
are from two distinct emission regions. 

The existence of two thermal components is also supported by a time-dependent spectral analysis 
of the X-ray data. The PL component derived from the BB+PL model shows an unreasonable 
hard-to-soft evolution feature, with the photon index continuously increasing from ~ 2 to 3.0 since 
t > 800 seconds (Fig. 2). This implies a very peculiar electron distribution. The PL photon index 
derived from the BB+BB+PL model, on the other hand, does not show a significant evolution, with 
an average index of 2.2 ± 0.4 (Fig.2). The unexpected softening derived in the BB+PL model is 
therefore an artifact due to the contribution of the soft photons in the high-T thermal component. 

Our time-resolved spectral analysis reveals the evolution of the two thermal components. Figure 
3 displays the temporal evolution of the temperatures and the effective blackbody radii of the two 
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components. The high-T component expands rapidly with Rh oc t and cools rapidly as kTn oc 
t~ - 87 , while the low-T component expands more slowly with R L oc t 38 and only cools mildly as 
kTi ~ t -0 - 05 . We note, however, that the peak flux of the low-T component is right at the low 
end of the XRT energy band where the absorption is very strong. It would be difficult to robustly 
measure the temperature evolution of this component. In Figure 4 we show both the unabsorbed 
and absorbed lightcurves of the two thermal components along with the non-thermal hard X-ray 
lightcurve in the 5-10 keV band. It is evident that the two components evolve independently. While 
the observed flux of the low-T component increases with time, the flux of the high-T component 
traces the non-thermal component. This suggests that the high-T component is directly related with 
the GRB jet. 

Numerical simulations of a relativistic jet propagating in the stellar envelope of a collapsing star 
reveals a hot cocoon surrounding the jet, 4 which is the waste heat from the jet enclosed by the 
dense envelope. When the jet, which may be mildly relativistic inside the star, first penetrates 
through the stellar envelope, the associated cocoon also breaks out and their thermal photons leak 
out. This hot thermal ring surrounding the jet then subsequently expands, both forwardly and 
in sideways, as the rest of jet bores its way out, The cocoon luminosity should roughly trace the 
non-thermal emission luminosity of the jet. This picture is in good agreement with the data of the 
high-T thermal component. Combined with the low-T component that has been interpreted as the 
shock breakout emission of the global supernova, 9 our result clearly suggests that there exist three 
components during the explosion of a dying massive star. Besides the more isotropically expanding 
supernova ejecta, there are indeed a collimated jet breaking through the envelope and producing the 
non-thermal emission and a hot thermal cocoon component associated with it. 

Assuming that the radius of the low-T component Rl ~ 3.2 x 10 12 cm is also roughly the radius of 
the exploding star, the high-T component then opens a very small solid angle of AQ ~ (R H /R L ) 2 ~ 
2.4 x 1CT 6 . The radius of the star may be smaller than R L , but in any case, it is unlikely that this 
solid angle reflects the solid angle of the relativistic jet, which should be very broad according to the 
late afterglow data 15,21 and population studies. 22,14 A possible picture is that the cocoon surrounding 
the jet forms a thin ring with a small solid angle. This is much narrower than the predicted values 
in numerical simulations. 4 However, XRF 060218 has a much lower luminosity than the modelled 
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Table 1. The comparison of the spectral fitting results with different models 
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typical collapsars. With much less waste heat surrounding an under-luminous jet, the cocoon could 
be in principle much thinner. 
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Figure 1 The observed spectra of the XRT data accumulated in the first orbit (150-2740 seconds 
since the burst trigger) and of the contributions from the two thermal components derived from the 
BB+BB+PL model fitting. 

Figure 2 The evolution of the power law photon indices (T) in the BB+PL model (open circle) 
and the BB+BB+PL model (solid circle), respectively. While the BB+PL case shows a strong 
evolution of T, no apparent evolution is revealed in the BB+BB+PL case, with an average T ~ 
2.2 + 0.4. 

Figure 3 The evolution of the blackbody emission radii (left) and temperatures (right) of the 
low-T (open circles) and high-T (solid circles) components, respectively, derived from the time- 
resolved spectral analysis. The solid and dashed lines are the regression lines for the high-T and 
low-T components, respectively. In our spectral analysis we first divide the XRT data into 50 
temporal segments. Each segment covers ~ 50 seconds. We generally require that the spectral bins 
are > 100 doe each segment. Some segments have a significantly smaller number of spectral bins. 
We then merge these segments with their adjacent ones. We finally obtain 32 temporal segments. 
The error bar in the time axis marks the time interval of each segment. We fit the spectrum with 
the BB+BB+PL model. In the time-resolved spectral analysis the N-r of the host galaxy is taken 
as N H = 0.64 x 10 22 cm -3 without considering its temporal evolution. 

Figure 4 The unabsorbed (a) and absorbed (b) light curves of the two thermal components 

derived from the time-resolved spectral analysis with the BB+BB+PL model and their comparisons 

with the light curve of the non-thermal component in 5-10 keV band (c). The fluxes are in units of 
iq-9 er g cm -2 g -i_ 
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Figure 2. 
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Figure 3. 



10 




Figure 4. 
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